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Comment on "Excitons in Molecular Aggregates 
with Levy Disorder: Anomalous Localization and 
Exchange Broadening of Optical Spectra" 

In their Letter [lj], Eisfeld et al. predicted exchange 
broadening and blue-shift of the absorption band as well 
as a nonuniversal disorder scaling of the localization 
length of excitons in J-aggregates with Levy disorder 
(a = h), which they contrasted with the previously ana- 
lyzed Gaussian and Lorentzian case 0] ■ The observations 
were explained by chain segmentation due to outliers, and 
its interplay with localization of states in effective poten- 
tial wells created by typical random site energies. We ar- 
gue that the previously known theory Q does not break 
down and thus anticipates the properties of absorption 
band investigated in the commented work. 



FIG. 1: a) JVioc scaling in a chain with Levy disorder. In- 
set: renormalized N\ oc distributions for different a values, b) 
Power law fits for Levy and Gaussian disorder; logarithmic 
scales should be used for correct results, c) Absorption band 
vs. bare disorder HWHM for Levy distribution. 

The nonuniversality of the iVi oc distribution shown in 
Fig. 3 of [l| does not arise from the properties of Levy 
disorder, but it is an artifact of considering a constant 
range of energies e € [—2.1, —1.9] J for comparison under 
different scale parameters a, instead of an interval ad- 
justing to disorder-induced scaling and shifts of absorp- 
tion band, as previously done for other disorder types 
(i.e. e € [—0.1, 0]) @. A cr-independent energy range im- 


plies that for low a values a too wide variety of states is 
included in the calculation. Consequently, the standard 
deviation 5N\ oc disproportionately increases and the ra- 
tio 5Ni oc / Ni oc is no longer constant. Using the correct 
scaling, e = ef, + (5e + 0.23)<7 2 / 3 , we have obtained univer- 
sal scaling of the Ni oc distribution with iVi oc ~ (J/ct) - 333 
and 5N loc ~ (J/ct) 324 (Fig.QJi). It breaks down only 
in the limit of small and large a due to the conver- 
gence towards N\ ac = 134 (unperturbed eigenstates) and 
N\ ac = 1, respectively. We have also recalculated the 
scaling of absorption band HWHM, particularly ensuring 
good accuracy for low a values, where the band becomes 
very narrow (Fig. lb). The obtained exponent 0.657, in 
contrast to 0.6 in the Letter, is in agreement with the the- 
ory in 0. We also note that in Fig. 2 of the HWHM 
should be compared not with a but with bare disorder 
HWHM, equal \J2 In la for Gaussian and 0.224cr for Levy 
distribution, in the latter case revealing the exchange 
broadening for all a values considered (Fig. lc) — not just 
[0,0.6] J. 

The above scaling relations result from the General- 


ized Central Limit Theorem, a* 


(1-q)/c 
loc 


and thus 


work for any a. The formula for the site-averaged disor- 
der strength carries the information about the type of its 
distribution (including heavy-tailedness for a < 2). The 
scaling N* ~ ( J / a) a ^ 1+a ^ results from the equilibrium 
between its typical value experienced by the absorption 
band states and the typical energetic cost of localizing 
them. From the universality of N\ oc distribution follows 
N\ oc ~ N*; together with the above relations, the univer- 
sality explains also the scaling HWHM ~ a{N*) {1 - a V a . 
On the other hand, segmentation (already present in 
Lorentzian disorder [2j) and localization in potential wells 
are just microscopic mechanisms realizing this statistical 
theory. The segmentation produces additional peaks in 
the absorption spectrum, predicted by the Letter, which 
for sufficiently low a can destroy the HWHM scaling al- 
ready in the intermediate disorder regime. 
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